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, Abstract 

r-j- We argue, that the multi-particle scattering amplitudes in N = 4 SUSY at large N c and in 

the multi-Regge kinematics for some physical regions have the high energy behavior appearing 
ON ' from the contribution of the Mandelstam cuts in the corresponding t-channel partial waves. The 

Mandelstam cuts correspond to gluon composite states in the adjoint representation of the gauge 
group SU(Nc). The hamiltonian for these states in the leading logarithmic approximation coin- 
cides with the local hamiltonian of an integrable open spin chain. We construct the corresponding 
f») ■ wave functions using the integrals of motion and the Baxter-Sklyanin approach. 

<D ' 

1 Introduction 

At high energies s ^> — t in QCD the elastic scattering amplitude for the process AB — ► A'B' in the 
■^j- ' leading logarithmic approximation (LLA) 

"sj" ■ a s Ins — 1 , a s < 1 (1) 



> 



, has the Regge form [T] 

o. 

O! A 2 ^ 2 = 2g5 XAXA ,TX A ,^^gT c BB ,5 XBXB ,, t = -f . (2) 



Here T c are the generators of the gauge group SU(N C ), X r are the particle helicities and j(t) = l+oj(t) 
■ is the gluon Regge trajectory for the space-time dimension D — 4 — 2e 

^ : »Hf) = I d^k *— « - a ( In ^ - i V (3) 



(2tt) 2 J k 2 {q- k) 2 V \ l 



e 



In the framework of the dimensional regularization the parameter fj, is the renormalization point for 
the 't Hooft coupling constant and 

( 4 ^" 7 ) e , 7 = "V>(1), (4) 

where 7 = — is the Euler constant and ip(x) — (lnr(x))'. The gluon trajectory j(t) was calculated 
also in the next-to-leading approximation in QCD [2] and in the SUSY gauge models [3]. 
In LLA gluons with momenta k r (r=l,..,n) are produced in the multi-Regge kinematics 

s = (p A +Pb) 2 > s r = (kr + fc r _i) 2 > -t r = ql, k r = q r+ i - q r , (5) 



*The talk given at the memorial Alcxei Zamolodchikov conference, June 21-23, 2008, Moscow, Russia. 
^Supported by the Maria Curie Award and grants RFBR 07-02-00902-a, RSGSS 5788.2006.2. 
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where the amplitude has the factorized form 



"(-4?) M-fi) s <"(-£+i) 
A 2 ^ 2+n = 2 S 5 XAXAI gTZ A ,^^gC^q 2 ,q 1 )e;(k 1 )T c % i ^^...^±l gT ^ S\\ B , . (6) 

<7l 12 Qn+l 

Here C IJi (q 2 ,qi) is the effective Reggeon-Reggeon-gluon vertex. In the case when the polarization 
vector e^ki) describes the gluon with a positive helicity in its cm. system with the particle A' one 
can obtain [4] 

C^C,(q 2 ,q 1 )e;(k 1 ) = V2^, (7) 

where the complex notation q = q x + iq y for the two-dimensional transverse vectors q was used. 

The elastic scattering amplitude with vacuum quantum numbers in the i-channel is calculated in 
terms of the production amplitude A 2 ^ 2+n with the use of the s-channel unitarity [1] . In this approach 
the Pomeron appears as a composite state of two Reggeized gluons. It is convenient to present the gluon 
transverse coordinates in the complex form together with their canonically conjugated momenta [HE] 

* ■ . d * . d 

Pk = x k + tyk , p k = x k - Wk , Pk = 1^— , p k = »a"Tr ■ ( 8 ) 

dp k dp* k 

In this case the homogeneous Balitsky-Fadin-Kuraev-Lipatov (BFKL) equation for the Pomeron wave 
function can be written as follows pQ 

E*(p L ,fa)=H i3 V(?L,fa), A = -^-^mmE, (9) 

where A is the Pomeron intercept entering in the asymptotic expression for the total cross-section 
a t ~ s A . The BFKL Hamiltonian has a rather simple operator representation [5] 

Ha = In \ PlP2 \ 2 + J^(k \ P i 2 \ 2 )piP* 2 + -^-(ln \ P i 2 \ 2 )p\p 2 - 4V(1) (10) 
PiP 2 P1P2 

with P12 = p\ — p 2 . The kinetic energy is proportional to the sum of two gluon Rcgge trajectories 
ui(~\pi\ 2 ) (i — 1,2). The potential energy ~ In \pi 2 \ 2 is obtained by the Fourier transformation from 
the product of two gluon production vertices C^. This Hamiltonian is invariant under the Mobius 
transformation [B] 

a Pk + b 

Pk~* —j, (11) 

cpk + a 

where a, b, c and d are complex parameters. The eigenvalues of the corresponding Casimir operators 
are expressed in terms of the conformal weights 

1 . n „ 1 . n , . 
m = — hi^H — , m = — h iv , (12 

2 2 2 2' y ' 

where v and n are respectively real and integer numbers for the principal series of unitary represen- 
tations of the Mobius group SL(2, C). The eigenvalues of Hi 2 depend on these parameters [B] 

E m ,fh = VK m ) + -m) + ip(m) + tp(l - m) - 4-0(1) . (13) 
The Pomeron intercept in LLA is positive 

A = 4— iV c ln2>0 (14) 

and therefore the Froissart bound at < c In 2 s for the total cross-section is violated |T] . To restore the 
broken s-channel unitarity one should take into account the contributions of diagrams corresponding 
to the ^-channel exchange of an arbitrary number of reggeized gluons in the i-channel. The wave 
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function of the colorless state constructed from n reggeized gluons can be obtained in LLA as a 
solution of the Bartels-Kwiecinski-Praszalowicz (BKP) equation [7] 

E^ = H^^f, A=-^^ min£ (15) 

47T 

In the N c — ► oo limit the color structure is simplified and the corresponding Hamiltonian has the 
property of the holomorphic separability [8 

n 

(°) = J2 H Kk+1 = ftW + h^> , [hP\ hP>] = . (16) 

k=l 

It is a consequence of the similar property for the pair BFKL hamiltonian H\2 (fT0|) and the energy 
The holomorphic Hamiltonian in the multi-color QCD can be written as follows (cf. (jTDJ) ) 

^ (0) = E h W > h x°2 = MP1P2) + ^ (lnpia)Pi + j 2 (lnpi 2 )p 2 - 20(1) , (17) 

where VK 2 -) = ( m r( x ))'- As a result, the wave function W has the holomorphic factorization [5] 

* = E pn ) ^(p*, , (18) 

which in the case of two-dimensional conformal field theories is a consequence of the infinite dimen- 
sional Virasoro group. Moreover, the holomorphic hamiltonian ft/ ) is invariant under the duality 
transformation [5] 

Pi -> Pm+i ~> Pi+i > (19) 

combined with its transposition. 

Further, there are integrals of motion q r commuting among themselves and with [51 110] : 

*r 0) = E Pk 1 k 2 Pk2k 3 --Pk 2 k 3 --Pk r k 1 Pk 1 Pk2---Pk r , [q r ,h]=0. (20) 

fcl<fe 2 <---<fer 

The integrability of the BFKL dynamics in LLA was established in Ref. [TO]. This remarkable 
property is related to the fact that h coincides with the local Hamiltonian of an integrable Heisenberg 
spin model [11]. Eigenvalues and eigenfunctions of this hamiltonian were constructed in Refs. [P21H3] 
in the framework of the Baxter-Sklyanin approach [14] . 

In the next-to-leading approximation the integral kernel for the BFKL equation was constructed 
in Refs. [21 HFJ- In QCD the eigenvalue of the kernel contains the Kroniker symbols S n> o and S n ^2 
but in N — 4 SUSY it is an analytic function of the conformal spin and having the property of the 
maximal transcendentality [31 116j . This extended supersymmetric theory appears in the framework 
of the AdS/CFT correspondence [T71 [TH] [19] . It is important, that the one-loop anomalous dimension 
for twist- 2 operators in N — 4 SUSY is proportional to the expression -0(1) — i'ii — 1): which is related 
to the integrability of evolution equations for the quasi-partonic operators in this model [20]. The 
integrability persists also for some operators in QCD [21]. The maximal transcendentality principle 
suggested in Ref. |16j gave a possibility to extract the universal anomalous dimension up to three loops 
in N = 4 SUSY [22] [23] from the corresponding QCD results [22 ]. The integrability of the TV = 4 
model was verified also for other operators, large coupling constants and in higher loops [25l [26l [27] , 
The asymptotic Bethe ansatz and integrability allowed to calculate the anomalous dimensions in four 
loops [2H]- The result is in an agreement with the next-to- leading BFKL predictions after taking into 
account the wrapping effects [33]. The maximal transcendentality was helpful for finding a closed 
integral equation for the cusp anomalous dimension in this model [301 131] with the use of the 4-loop 
result [32] . 
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There is another region of investigation, in which remarkable properties of the N=4 SUSY are also 
found. Namely, Bern, Dixon and Smirnov (BDS) suggested a simple ansatz for the gluon scattering 
amplitudes in this model 33J . This ansatz was verified for the elastic amplitude in the strong coupling 
regime using the AdS/CFT correspondence [34]. But the BDS hypothesis does not agree in this regime 
with the calculation of the multi-particle amplitude [35j . The property of the conformal invariance 
of the BDS amplitudes in the momentum space was discussed in Ref. [36] and the relation with 
the Wilson loop approach was suggested in Ref. [37] generalizing the results of the strong coupling 
calculations of Ref. 34 . The BDS amplitudes A n for n > 6 in the multi-Regge kinematics do not 
have correct analytic properties compatible with the Steinman relations [38]. It is a consequence of 
the fact, that these amplitudes do not include the Mandelstam cuts [38 . This cut contribution was 
obtained from the BFKL-like equation for the amplitude with the i-channel exchange in the adjoint 
representation of the gauge group [38] . This equation was solved in LL A and the two-loop expression 
for the 6-point scattering amplitude in the multi-Regge kinematics was derived [39j . Recently the 
two-loop correction was calculated numerically for some values of external momenta in an agreement 
with expectations based on the Wilson loop approach 40J. 

In this paper we demonstrate, that in the multi-color limit for the production amplitudes the 
contributions of the Mandelstam cuts generated by the multi-Reggeon i-channel exchange can be 
expressed in terms of the solution of the BKP-like equation for the composite states of several reggeized 
gluons in the adjoint representation. It turns out, that in LLA the corresponding Hamiltonian coincides 
with the local Hamiltonian of an integrable open Heisenberg spin chain. These results partly were 
presented at the conferences [JT] H2] . 

2 Mandelstam cuts 

A planar amplitude for the production of two gluons in the multi-Regge kinematics s ^> \s%\ ~ | S2 1 ~ 
l s 3| 3* |ti| ~ 1^1 ~ 1^3 1 has the multi-Regge form almost in all physical kinematical regions. But in the 
physical region where Si, S3 < 0; s > 0, S2 > the amplitude contains also the Mandelstam cut [43] 
in the angular momentum plane j'2 of the crossing channel ti = —q 2 [38] in the adjoint representation 
of the color group. The cut appears as a result of the exchange of two reggeized gluons with the 
momenta p\ — k and P2 = q — k, respectively [39J (see Appendix A for more details). In the region 
si, S3 < 0; S2 > the integrals over the Sudakov variables a = kpA/pAPB and (3 = kps/pAPB do not 
vanish as in other regions because the integrand contains singularities situated above and below the 
corresponding integration contours in an accordance with the Mandelstam requirements [43]- These 
singularities lead to simultaneous discontinuities of the amplitude in the invariants S2 and s. 

For the planar amplitude with six external particles only diagrams with two reggeons in the ti- 
channel give a non-vanishing contribution because for a larger number of reggeons the Mandelstam 
conditions for singularities in other Sudakov variables are not fulfilled. However in the case of a 
larger number of external particles the exchange of several reggeons with momenta pi gives also a 
non-vanishing contribution to the amplitude constructed from planar diagrams. For the Mandelstam 
cut resulting from an exchange of n reggeons one needs at least k=2+2n external particles to have 
simultaneous singularities in upper and lower complex semi-planes for the Sudakov parameters a[, (3[ 
of the reggeon momenta pi, as it is demonstrated in Appendix A. 

Let us discuss such composite state of n reggeized gluons in the adjoint representation (cf. a similar 
approach for the simple case n = 4 in Ref. [39] ) . One can write the homogeneous BKP equation for 
its wave function described by an amplitude with amputated propagators in the form (see Appendix 
A) 

H^ = EV, A n = -^E. (21) 
Here H is a redefined hamiltonian obtained after subtraction of the gluon Regge trajectory ui(t) 
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containing infrared divergencies. Namely, the Regge trajectory of the composite state is [41] [42] 

Un(t) =a(j~ ln Jpj + A » . a = (4^e-T) £ , (22) 

where A„ is the infrared stable quantity expressed in terms of the energy E. 

The hamiltonian H in the multi-color limit can be written in the holomorphically separable form 
(see Appendix A) (cf. [3S]) 

n— 1 71 

H = h + h*, h = In^|2 + J2 h r,r+l , 9 = , (23) 

^ r=l 1 

where the pair hamiltonian h rr+1 is transposed to the corresponding unamputated operator (|17p 

h\ r+1 = ln(p r p r+ i) + p r ln(p r r+ i) — + p r+i ln(p r r+1 ) — h 27 . (24) 

Pr ' Pr+1 



It is seen from eq. (|23[) that the holomorphic hamiltonian for the composite state in the adjoint repre- 
sentation differs from the corresponding expression for the singlet case (fT7| after its transposition 
only by the substitution 

h^-ln^, (25) 
q 2 

which is related to the fact, that the planar Feynman diagrams have the topology of a strip and the 
infrared divergencies in the Regge trajectories of the particles 1 and n are not compensated by the 
contribution from the pair potential energy V n ,\- 

It turns out, that the eigenvalues E do not depend on \q\ 2 due to the scale invariance of H, as it 
will be demonstrated below. As a result, the ^-dependence of u> n {t) is the same as in the gluon Regge 
trajectory. 

The transposed holomorphic hamiltonian is related to the initial hamiltonian by the similarity 
transformation 

h* = {f[P^j h (j[ Pr J ' (26) 
which leads to the following hermicity property of the total hamiltonian H 

h+ = (n h 2 ) h (n \ Pr \^ . (27) 

The last relation is compatible with the normalization condition for the wave function 

/n— 1 n n 

l[cFp r **l[\ Ps \- 2 *, Y d P. = q. (28) 
r=l s=l s=l 

Using the duality transformation (cf. [9]) 

Q 

Pi = ^0,1 , Pr = Zr-l,r , 9 = Z0,n , Rr,r+1 = = id r , (29) 

the holomorphic hamiltonian can be rewritten as follows 



n-l 

h = In + Y: K, r+1 , (30) 

: n r= i 
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where 

ft*,r+l = 2 ln ( 5 r) + T. + T. + + 1 Z r -l,r) + ^ . (31) 

Here and later we neglect the pure imaginary contribution 2 ln(i) because it is cancelled in the total 
hamiltonian H. Note, that for the colorless composite state and q — zo^ n — the transformation (|29|) 
is indeed reduced to the usual duality substitution of Ref. [9] . 
To simplify h one can use the relations 19] 

In 9 = - Inx + ^(ip(xd + 1) + tpi-xd)) , 

\n(x 2 d) =lnx + -{^{xd) + ip(-xd + 1)) , 

In d = ln(x 2 d) - 2 In x + 4 - . (32) 

o x 

Then h l T r+1 can be presented as follows 

htr,r+l = l n ( z r,r+l^r) + \xy[z r _i r d r ) — ln^r,r+l — ln2 r — l,r + 2 7 . (33) 

Further, by regrouping its terms we can write the holomorphic hamiltonian in another form 

h = -21nz ,„ + ln(zo,i<9i) + ln(4-i,„9„-i) + 2 7 + ^ h' r r+1 , (34) 

r=l 

where 

K,r+1 = M Z r,r+l d r) + M Z r,r+l d r+l) ~ 2\nZ^ r+1 + 2j 

= ln(d r ) + ln(5 r+ i) + — In z r , r+1 d r + — !— In z r , r+ i d r+ i + 2j . (35) 

Or Or+1 

The pair hamiltonian h' rr+l coincides in fact after the substitution z r — > p r with the corresponding 
hamiltonian in the coordinate representation (|17p acting on the wave function with non-amputated 
propagators. 

In particular, for n = 2 one obtains (cf. [55] ) 

h= -21n^o,2+M2o,i9i)+ln(42^i) + 2 7- (36) 
It is important, that h (j3"4"|) is invariant under the Mobius transformations 

az k + b 

z k -> —7 (37) 

czfc + a 

and does not contain the derivatives do and 9 n . Therefore we can put 

z — , z„ = oo , (38) 
which leads to the simplified expression for h 

n-2 

h^ti = lnizfdi) + ln(3„_i) + 2j + J2 K,r+i • (39) 

r=l 

To return to initial variables in the final expression for the wave function one should perform the 
following substitution of Zfc 

^fPf^^N^-- (4o) 
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According to the above representation (|30|) for h, its transposed part h* can be obtained from h 
by the similarity transformation which can be written in terms of h! as follows 

(n-2 \ _1 {n-1 \ 

n Zr ' r + i ) h> zi (n Zr - r + i J ' ( 4i ) 

which is compatible with the following normalization condition for the wave function in the full two- 
dimensional space 

V^IIr^w 2 - (42) 



On the other hand, from the expression (I39p for h! we obtain another relation for h! 



-l 



h>t = ( n>j h> [J\ dr ) ' (43) 

corresponding to the second normalization condition for \l/ compatible with the hermicity properties 
of the total hamiltonian 

/n — 1 n — 1 

n^^n^i 2 *- ( 44 ) 
r=l r=l 

By comparing two above relations between h! and h n one can conclude (cf. [10]), that the operator 

n—2 n-1 

A' = J] z s , s+ i [| o> r (45) 

s — l r—1 

commutes with the holomorphic hamiltonian 

[A',h'] = 0. (46) 

3 Integrable open spin chain 

Let us verify, that the holomorphic hamiltonian h' (|39[) also commutes with the differential operator 
D(u) being the matrix element T22 of the monodromy matrix (cf. [10] 'l 

T i u ) = ( n [ ul nul ) = L 1 (u)L 2 (u)...L n . 1 (u) , (47) 



C(u) D{u 

where the L-operator is defined by the relation 



L r (u)=[ U + 'lf ) . N x, 



-iz^d r u — iz r d r 

To prove the commutativity of h! and D(u) one can use the following relation 

[L k (u) L k+1 (u), h' kk+1 ] = -i (L k (u) - L k+1 (uj) , (49) 
valid due to the Mobius symmetry of the pair hamiltonian 

[M k! k +1 ,h' k>k+1 \ = , M k , k+1 = M fc + M fc+ i (50) 
and the commutation relation (see [9]) 

[K.k+n [Ml k+1 ,N kik+1 ]] = 4A\ fc+1 , N k<k+1 = M k - M fe +i ■ (51) 
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The last relation is a consequence of the fact, that the operator Nk,k+i has non-vanishing matrix 
elements only between the states \m k ,k+i > and \m,k,k+i ± 1 > in the representation, where the 
Casimir operator of the Mobius group is diagonal 

M l k+i\ m k,k+i >= m kik+1 (m kik+1 - l)|mfc,fe + i > . (52) 



In this representation the commutation relation (|51| is reduced to the recurrent relation for the 
eigenvalues e(m k k +i) of the hamiltonian h' k k _^_i 

e(m + 1) — e(m) = 2/m , (53) 

fulfilled due to the well known representation of e(m) 

e(m) = ip(m) +ip(l ~ m ) + 2 7' ( 54 ) 
Relation P9"| leads to the equality 

n-2 

[T(u), ^ K,r+i\ = iL2(u)L 3 (u)...L n _i(u) - iLi(u)L 2 (u)...L n - 2 (u) . (55) 

r=l 

On the other hand, one can easily verify, that 

[T 22 (u)M49i) +lnd n -i] = (0, 1) [T(i0,]n(*?&) +ln$ l _i] ( J ) 

= -i (0 , 1) (L a (u)i3(«)...£»-i(«) - Li(«)La(«)...Ln-a(«)) ( J ) - ( 56 ) 

which proves that the differential operator D{u) = T 22 {u) is an integral of motion 

[D(u),ti] = 0. (57) 

Thus, our hamiltonian is the local hamiltonian for an open integrable Heisenberg spin model with the 
spins which are generators of the Mobius groupQ- 

With the use of the following decomposition of the L-operators 



L r (u)=^ Q u j + ^_ zr j(z r ,l)id r (58) 
one can construct the matrix element T22 — D(u) in an explicit way 

n-l 

D(u) = Y / n n - 1 - k q' k , (59) 



fe=0 

where 

q'o = 1 , q[ = ~i ^2 z r d r . (60) 



n-l 



In a general case the integrals of motion q' k are given below 

fc-1 k 



q'k = - Zr ^ II Z r*,r 3+1 Y[id rt ■ (61) 

0<ri<r 2 <...<r fc <n s=l t=l 



1 l thank L. D. Faddeev for the fruitful discussion in which he suggested, that the operator D(u) could be an integral 
of motion for this open spin chain. 
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In particular, we obtain, that q f n _i is proportional to the integral of motion A! (145ft 

n—2 n—1 

<&_i - -i n ~ x z x J] z s , s+1 l[d t = -i n - x A' . (62) 

s=l t=l 

Note, that one can parameterize the monodromy matrix in another form 

Tin) = ( j ° (u } + ? (u) , ,. A ) , i ± («) = *(«) ± «,(«) • (63) 



jo(u)-j 3 (u) 

In this case the Yang-Baxter equations for the currents have the Lorentz-invariant representation 9 

[?», = j^-r {f{u)f{v)-j»{v)f{u)) . (64) 

2{u — v) 

Here e^pa is the antisymmetric tensor in the four-dimensional Minkowski space and £1230 = 1 , g^v = 
(1,-1,-1,-1). 

In particular, we obtain from the Yang-Baxter equations the relation 

bo(w) - j3{u),j {v) - j 3 (v)} = \jo(v),h(u)] - \jo(u),j s (v)] = (65) 
and therefore the integrals of motion q' k are independent operators and commute each with others 

[q' k ,ql\=0. (66) 

4 Composite states of two and three gluons 

In the case n — 2 we have only one non-trivial integral of motion 

q[ = -iz\ di . (67) 
Taking into account the normalization condition for the eigenfunction in the two-dimensional space 

\\m 2 = /^i*! 2 ' (68) 

we find the orthonormalizcd and complete set of eigenfunctions 

H>W~=z 1 * {zl) * +m , m= —^-+ tv , m= —^--tu, (69) 

satisfying the single- valuedness requirement. Note, that using the substitution (I40[) one can reproduce 
the wave functions of two gluon composite states in the momentum space (sec 39 ) . 
For the case n = 3 the operator D(u) is given below 

D 3 (u) = u 2 - iu(zi di + z 2 d 2 ) + z\Z\ t i d\d 2 ■ (70) 

With taking into account the normalization condition 

d 2 Z\ d 2 z 2 
\zi\ 2 \z h2 \ 2 



*H = / l„ 121, 12 1^1 ■ 



one can search the holomorphic eigenfunction of this operator in the form 

4+m , ( z 2 



(72) 



The function f(x) satisfies the equation 

x(l - x)d 2 + (\ + m)(l - x)d + x] f = 0, x = - , (73) 
2 J Z\ 

where A is the eigenvalue of the operator z\Z\^ 2 d\d 2 . Two independent solutions of this equation can 
be expressed in terms of the hypergeometric function F 

h{x) = F(a u a 2 ; 1 + ai + a 2 ; x) , f 2 (x) = x ai+a2 F(-a 2 , -ai; 1 - ai - a 2 ; x) , (74) 

where the parameters a\ and a 2 are obtained from the set of equations 

ai + a 2 = + m , aia 2 = -A . (75) 

The solutions near the point x = 1 can be also expressed in terms of hypergeometric functions and 
are expanded as follow 

rn+i^+lV 1 ^!-* = r~r -( x ~ l ) (Mi -x)- V(i) m + ^(l + 0l ) + V(i + « 2 )) 

1 (1 + ai + a 2 J aia 2 

and 

r(- 0l )r(- 02 ) = J_ _ {x _ 1} (ln(1 _ x) _ _ ^ (2) + ^ (1 _ ai) + ^ _ aa)) 

1 (1 — ai — a 2 ) aia 2 
Analogously one can find the large- x behavior of the functions /i and f 2 

, _ r(a 1 +a 2 + l)r(a 2 -a 1 ) / . T( ai + a 2 + l)r(ai - a 2 ) 2 

M j| — " r(a 2 )r(i + a 2 ) 1 ^ + r(ai)r(i + 0l ) 1 j ' 

, _ r(l-a 1 -a 2 )r(a 1 -a 2 ) 2 r(l - Oi - a 2 )r(a 2 - fll ) oi 

/2(xj i— - r(-a 2 )r(i-a 2 ) ( x) + r(- ai )r(i- ai ) ( x > ' 

To construct the wave function with the property of the single- valuedness in the two-dimensional 
subspaces z\ and x we should write a bilinear combination of the functions fi (x) and the corresponding 
functions in the anti-holomorphic subspace fi(x*) taking into account that in the second pair of 
functions one should perform the substitution 

~ ~ ~ 1 — n /m~\ 

ai — ► ai , a 2 ^> a 2 , m — ► m = — h iv . (76) 

Due to the single- valuedness of the wave function near i = 0we obtain for it the following expression 

*=N 2w (^V , ^ m ,fh{x)=h{x)h{x*) + Cf 2 {x)J 2 {x*), (77) 



where the constant C should be fixed from the requirement, that the analytic continuation of <J/ in the 
neighborhood of the points x = 1 and x = oo leads also to a single- valued expression. The condition, 
that near x = 1 the terms proportional to \x — 1| 2 ln(l — x) ln(l — x*) are absent, gives the relation 

r(ai + a 2 + 1) T(aj +a 2 + l) T(l - oi - q 2 ) r(l - ai - a 2 ) _ 

T( ai )T(a 2 ) r(oi)r(a 2 ) r(-oi) r(-oa) r(-ai) r(-a 2 ) ~ ■ 1 j 

Providing, that the constant C is fixed by this equality, the behavior of the total wave function at 
x — > 1 is simplified 

lim tf m>fa (£) ~ (V>(1 + ai) + i>(l + a 2 ) - V>(-ai) - ^(-o 2 )) |1 — ln(l - a;*) 

x— ► ! 
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+ (ip(l + Si) + + a 2 ) - V'(-Si) - ip{-a 2 )) |1 - x| 2 ln(l - i) . (79) 
Thus, the single- valuedness condition at x — > 1 leads to the additional equation 

cot(7rai) + cot(-7ra2) = cot(7rai) + cot(7ra2) • (80) 

A stronger constraint can be obtained from the single- valuedness condition for ^ at x —> oo. 
Indeed, its consequence for the bilinear combinations 

(-x)-°(-s*) _sr , (-x)- a2 {-x*)-™ 

leads to the relations 

ai - Si = iV 01 , &2 - S2 = N a2 , (81) 
where N ai , N a2 are integers. Further, the absence of the interference terms 

(-x)- ai (-x*)-^ , {-xY a2 (-x*)-^ 



is fulfilled due to the above relation (|77|) for C. 

One can write the integral representation for the wave function satisfying the above constraints 



\y\ 2 ' \y-%j \y -x* j z x 

where the integration is performed over the two-dimensional plane y. Note, that the integrand has no 
ambiguity in the points y = 0, l,x due to the derived relations between ai,Si and 02,02- Moreover, 
the function ^> near the points x — 0, 1, 00 can be presented in terms of the sum of products of above 
hypergeometric functions. 

There is another basis for the holomorphic solutions 

z 2 ) = zl 1 z% 2 F(oi, -a 2 , 1 + 01 - a 2 ; — ) , 

z 2 

^ 2 (zi,z 2 ) = z\ 2 z^ 1 F(a 2: -at, 1 + a 2 - 01; — ) (83) 

z 2 

allowing to construct an equivalent representation for the total wave function 'J. Note, that these 
functions can be written in terms of the Mellin-Barnes integrals 

J -too r ( a i -a 2 + l + s) 

M*>»)~r T{a ^; )T{ -\ + A s }\ { - s) (-z ir ^(-z 2 r- ds . (84) 

J -too r(a 2 -ai + l + s) 

Here it is assumed, that the poles of r(— s) are situated to the right from the integration contour 
whereas all other poles lie to the left of it. 

5 Hamiltonian and integrals of motion 

The holomorphic hamiltonian for composite states of two reggeized gluons can be written as follows 
h = Hzfdi) + ln(<9i) + 2 7 = ^(zi<9i) + ip(-Z!di) + 2 7 . (85) 
Acting by h on the function z\ we obtain 

hz{ = e{8) z{ , e(5) = $(8) + ${-8) + 2 7 . (86) 
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In the case of wave function ([69)1 satisfying the single- valuedness and orthonormality conditions in the 
two-dimensional space one derives the following expression for the total energy |39j 

E m ,m =£ m + (fn , £m = V>( _ 2 + m ) + ^(^ ~ m ) + 2 7 • ( 87 ) 

Note, that it does not coincide with the corresponding result (fT3")l for the Pomeron state. 
The holomorphic hamiltonian for composite states of three gluons has the form 

ti = ln(z^i) + ln(9 2 ) +ln(z? i2 9i) + In^ft) - 21nz 1 , 2 + 4 7 . (88) 

In the region 

Z\ < z 2 (89) 
it is a sum of two independent pair hamiltonians 

ti = ip(zidi) + ip{-zxdx) + ij{z 2 di) + ip{-z 2 di) + 4 7 . (90) 

Because the limit x — z 2 j z\ — ► oo in solution (|T2[) corresponds to this kinematics, we obtain 

e = e(oi) + e(oa) , (91) 

where ai and a 2 are parameters of the three-gluon composite state (see (|75)l ). The eigenvalues of the 
integrals of motion are also expressed in terms of these parameters. Due to the normalizability condi- 
tion these quantities together with the parameters oi, 02 of the wave function in the anti- holomorphic 
space should be chosen as follows 

ai = w ai + — , a 2 = iv a2 + — , 

ai = w a x - -y- 1 02 = w a . 2 - ~y , (92) 

where v r are real and n r are integer numbers. 
Note, that 

v = v ai +Va 2 , n = n ai + n a2 , aia 2 = -A , a x a 2 = -A (93) 
and the eigenvalues of two integrals of motion q' k can be obtained as coefficients of the polynomials 

P2(u) = (u — iai){u - ia 2 ) , ^2(1*) = (u — ia L ){u - ia 2 ) ■ (94) 

Generally for the composite state of n reggeized gluons the situation is similar. Namely, the 
holomorphic wave function in the region 

zi < z 2 < z 3 < ... < z„_i . (95) 

is factorized 

n-l 

*o ll( ^ l ... 1 . B _ 1 = n < r ■ ( 9e ) 

r=l 

The energy for this solution is the sum of the particle energies 

n-l 

e = J2<«r). (97) 

r=l 

The eigenvalues of integrals of motion q' k can be expressed in terms of the coefficients of the 
polynomial 

n-l 

Pn(u) = ]]_(u-ia r ). (98) 

r=l 
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Due to the condition of the normalizability the parameters should have the form 

a r = iv T + , (99) 

where v r is real and n r is an integer number. The energies and eigenvalues of the integrals of motion 
in the anti-holomorphic space are given by the same expressions with the corresponding substitution 
of parameters 

a r -^a r = iv r — ^ . (100) 

The holomorphic wave function satisfies a set of differential equations following from the eigenvalue 
equation for the operator D(u) 

n-l 

D(u) *ai,o 2 ,...,o n _i = II ( M ~ iCl ^ *oi,o 2 ,...,o„_i • (101) 
r=l 

This equation can be solved with the use of the Taylor expansion 

= II ^ E (?) -El 7^) ( 102 ) 

r=\ s 2 =0 y ' s n -l 

where the coefficients c(s2, are calculated in a recurrent way. The recurrent relations obtained 

from the eigenvalue equations for different operators q' r are compatible due to their commutativity. The 
obtained solution has the singularities at Zki = 0. But, if we consider (n — 1)! functions * ail ,.. 0i 
obtained by all possible permutations of parameters a r and multiply them on the corresponding 
functions in the anti-holomorphic subspace, it is possible to construct the wave function having the 
single-valuedness property in two-dimensional spaces z r 

tf(2i,...,*„_l) = Cf {ii,...,i»-i}*a (ll a i!1 ,...,o 4n _ 1 *o 4l ,o, 2 ,...,o, n _ 1 ■ (103) 

{ii,i 2 ,...,i„-i} 

For this purpose one should adjust the coefficients Cr^ . ^ j in an appropriate way presumably 
without additional constraints on the parameters a r and a r . The composite state of n — 1 gluons has 
the following total energy 

n— 1 n— 1 

£ = e + e, e = ^ e(o r ) , e = ^ ■ (104) 

r— 1 r— 1 



6 Baxter-Sklyanin approach 

To find a solution of the Yang-Baxter equation for the open spin chain one can use the Bethe ansatz. 
For this purpose it is convenient to work in the transposed representation for the monodromy matrix 

r'M - ( fl(B f « ) - £1 (OUW.-H., (.) , (106) 

where the L-operator can be chosen as follows 

4(«)= ( U+ .i dr * r l % ). (106) 

r ' \ —lO r Zr U — lO r Z r J v ' 

The pseudo-vacuum state is defined as a solution of the equation 

£(«)*o = 0. (107) 
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It can be written in the form [11] 



n-1 



*o = [I z- 2 . (108) 



r=l 

2 



Note, that the function does not belong to the principal series of the unitary representations. 

As a result, the states constructed in the framework of the Bcthe ansatz by applying the product of 
the operators j\{u r ) to 

k 

*fc = nj+(«r)*0 (109) 

r=l 

are non-physical. Nevertheless, these states are eigenfunctions of the integral of motion 

= (jl{u) - jt(u)M = A(u)¥ k (110) 

providing that 

is a polynomial, which leads to a quantization condition for the Bethe roots Ut- If we parameterize 
this polynomial as follows 

n-1 

A(u) = Y[(u - leu) , (112) 
i=i 

the above defined Baxter function Q(u) can be calculated 

Here for generality we included the factor <p(u) which is an arbitrary periodic function 

<t>{u) = <t>{u + i) . (114) 

In the case of a finite number of the multipliers j\.{u r ) in the Bethe ansatz for the wave function 
the expression Q(u) is also a polynomial 

fe 

Q(u) = J[(U - U r ) . (115) 

r=l 

For such solutions the parameters ai = — ki — 1 are negative integer numbers satisfying the condition 

n-1 

J2h=k. (116) 

i=i 

The corresponding Baxter functions can be written as follows 

n— 1 ki maxt kt 

Q(u)=l[l[(u + tt)= J] (u + ipY P , (117) 
i=i t=i P =i 

where r p is the number of k t satisfying the condition k t > p. 

As it was mentioned above, the polynomial solutions for Q(u) are non-physical, because the cor- 
responding wave functions \& do not belong to the principal series of unitary representations of the 
Mobius group. We should find a set of non-polynomial solutions Q s (u) satisfying this physical re- 
quirement. 
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According to E. Sklyanin [14] the correct variables in which the dynamics of the Heisenberg spin 
model is drastically simplified are the zeroes b r of the operator B(u) = j^{u) entering in the mon- 
odromy matrix 

n — 2 n — 1 

B(u) = P n _a H(u- b r ) , P n _i = i d r , (118) 

fe=l r=l 

where the operators 6 r and P n ~\ commute each with others 

[b r ,b s ] = [b r ,P n - 1 ]=0. (119) 

It is convenient to pass from the coordinate representation z to the Baxtcr-Sklyanin representa- 
tion [12], in which the currents j+(w) and (j+(u))* (together with the operators b r , b* and P n -i, Pn-i) 
are diagonal. We denote the eigenvalues of the Sklyanin operators by b r ,b*. The kernel of the uni- 
tary transformation to the Baxter-Sklyanin representation is known explicitly for the cases n = 2, 
n = 3 and n = 4 [12] . For general n this integral operator can be presented as a multi-dimensional 
integral [TB"] . 

In the Baxter-Sklyanin representation the wave function in the holomorphic subspace can be 
expressed as a product of the pseudo-vacuum state in this representation Wo(P n -ii bx, 62, frn-2) 
and the Baxter functions Q(ut) 

n-2 

**(P n -l5 h, 6»_ 2 ) = f„"-T _m II *0(Pn-l, bl, fc«-2) , (120) 

k=l 

where the power-like behavior in the variable P n -i is in an agreement with the normalization condition. 
The analogous representation is valid for the total wave function 

n-2 

*+(£„_!;&!, ...,b n - 2 ) = P~^f~ m {P*_ x )-^-™ Y[Q(b r )^o(Pn-i;h,...,b n . 2 ) (121) 

fc=l 

with the use of the generalized Baxter function Q(u) being a bilinear combination of the usual Baxter 
functions in the holomorphic and anti-holomorphic subspaces 

Q{u)=J2<tQ s {u)Q t (u*). (122) 

s,t 

Here Q s (u) are different solutions of the Baxter equation with the same eigenvalue A(u). The coeffi- 
cients d Sy t are chosen from the requirement, that the function Q(u) satisfies the normalization condition 
everywhere including the points where the functions Q s {u) and Qt(u*) have the poles [T2l [T3] . For 
the periodic spin chain this condition leads to the quantization of the eigenvalue of the operator 
A(u) + B(u) although a simpler method of quantization is based on the requirement, that all Baxter 
functions corresponding to the same eigenvalue should have the same holomorphic energy [12] . In the 
case of the open Heisenberg spin model the situation is simpler and will be discussed below. 



7 Baxter-Sklyanin representation for two and three gluon states 

Let us consider the composite states constructed from two and three reggeons in the framework of 
the Baxter-Sklyanin approach. In the case n = 2 we have the following integral of motion in the 
transposed space 

D t (u)=j -j 3 = u-id 1 z 1 (123) 
and its eigenstates in accordance with the Sklyanin approach are given by the expression 

* 2 z x ~ H ■ (124) 
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The corresponding transposed hamiltonian is presented below 

/i* = ln(di^)+lnd 2 + 27. (125) 
Its eigenvalue calculated on the above cigenfunction 1 3>* is 

e ro = V(-i+m)+V(^-m) + 2 7 . (126) 

For the states composed from three reggeized gluons the transposed integral of motion in the 
holomorphic subspace is 

D\{u) = u 2 - iu{d 1 z x + d 2 z 2 ) + 8^2 Zl z h2 (127) 
and the operator j l + is given below 

3 \ = iu{d l + d 2 ) - d 1 d 2 z 12 = i{d! + d 2 ) (u - k) , (128) 

where 

^-^TTTir^. (129) 
di + d 2 

The operator j\ is easily diagonalized after a transition to the momentum representation, where 

id\ f Pl , P2 — Pi Ipi,p2 > *<9 2 /pi,p2 — P2 fpi,p 2 ■ (130) 
In this case the eigenvalue equation for jt_ has the form 

d d \ 

u(pi+p 2 ) - ipip 2 (-^ - Q^)j f = (Pi +P2)(u - bi)f, (131) 
where b\ is the eigenvalue of b\ . Its solution is given below 

\ —ih\ 

Pi 



f = x(pi+P2M)\j 2 J , (132) 

where \ is an arbitrary function of p\ + p 2 and bi . The dependence of \t* from p\ + p 2 is fixed by the 
normalization condition 

** ~ (pi +P2y ai ~ a2 . (133) 

On the other hand, the eigenvalue equation for the integral of motion in the momentum space can 
be written in the form 

PlP2 Wi ~ w) * (pi ' P2) = aia2 * (pi ' P2) • (134) 

Using the anzatz 

*(pi,P2) = (pi + P 2)- ai - a2 v(y) , y = — , (135) 

pi 

we obtain the following equation for the function rj{y) 

(y 2 d 2 + (ai + a 2 + 1) y d - ai a 2 ) V (y) = (-y 3 3 2 - 2y 2 8) V (y) . (136) 
There are two independent solutions of this equation 



£i r(fc + i)r(fc)r(i- ai )r(i-a 2 ) 
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-F(l-Oi,l-a 2 ,2;--) = ; rF(-0!,l- 01,1 + 02-0x5-2/) 

y 2/ r(i - a 2 ) r(i + ax) 

r '" 2 " F(-o 2 ,l-a 2 ,l + ai-a 2 ;-y) (137) 



r(i- ai )r(i + a 2 ) 



and 



1 ^ r(fc-ai)r(fc-a 2 )(-i) fc - 1 2/ - fc „ 

m{y) = + 2, Wi i TvrvtTvrTi vFri V (my + c k {a 1 ,a 2 )) 

ai a 2 f— T(k + l)r(fc)r(l - oi)r(l - a 2 ) 



fe=i 

= - ru ( y F _ai ' 1 - oi. 1 + as - ai; -j/) , 138 

1 (1 + a 2 - ai) 

where 

Cfc(ai, as) = ^(fc) + ip(k + 1) - - m) - ^(1 - fc + o 2 ) . (139) 

One can construct the bilinear combination of these solutions having the single-valuedness property 
at y = oo 

v(v) ~ vi(y)m(y*) + V2{y)m(y*) + c vi(y)m(y*) ■ (140) 

On the other hand let us use the above expression for ??i(y) and 772(2/) expressed in terms of the 
hypergeometric function regular at y = 0. To cancel the interference terms violating the single- 
valuedness condition at y — > in the above bilinear combination for rj(y) we should fix C as follows 

~ _ sin(ai7r) sin(a 2 7r) _ sin(ai7r) sin(a 2 7r) (141) 
7rsin((ai — a 2 )7r) 7rsin((ax — a 2 )-7r) 

Finally with the use of the integral representation for the hypergeometric function the wave function 
^P* in the momentum space can be written as follows 

00 = (pi + P 2)- ai - a2 (p*i +p^)- 31 - 32 m , (W2) 

where (f>(y) is given below 



«^)=/a^(- + l) ^ + 1) (l-t)" 3 -^!-**)" 3-1 (143) 



and satisfies the single valuedness condition in the y-space due to the quantization conditions (|92p . 

The transition to the Baxter-Sklyanin representation (it, u) corresponds to the Mellin-type trans- 
formation of <p(y) 

4>{u, 5) = / S y" m (y*r a <t>(y) = J <* 2 t (i - 1)"- 1 (i - rf^ 1 x (f) , (144) 

where 

- zit = zz/„ + — , -iu = iv u — . (145) 

Here v u is a real number and N u = 0, ±1, ±2, .... The function \ is given below 



\rt= lw y ~ iu(y * ri *(h + 1 ) 1 G^ +l) =/ "' ,/ ,l4(iJ 



The corresponding integrals can be calculated explicitly 

r(l + ox) T(iu)T(-iu- ai ) 



C\ = 7T ■ 



r(-ox) r(l -iu)T(l + iu + ai) ' 
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c 2 = I (i - tr (i - t*r t iu (t*r = tt 

Therefore we obtain for </>(u, u) the following expression 



r( fl2 ) r(i + m)r(-w-a 2 ) 

r(l - a 2 ) T(-iu) T(l + iu + 02) 



_ 7r 2 r(l + ai)r(a 2 ) r(m)r(i + iu) r(-m - 01) r(-m - a 2 ) 

^ U ' UJ r(-ai)r(i-o 2 ) r(-iu) r(i - iu) r(i + m + Si)r(i + m + s 2 ) ' 1 ' 

The inverse transformation corresponds to the Baxter-Sklyanin representation for the wave function 

// \ — iu / * \ — iu 
d 2 u0(u,u) l^-j (j|J , (148) 



where 



iu = iv u + , -iu = iv u , / d 2 u = / oV„ . (149) 

2 2 ^ • J '- 00 AT„=-oo 



One can interpret the wave function <f>(u, u) in the Baxtcr-Sklyanin representation as a product of 
the pseudo- vacuum state u u and the total Baxter function 

4>(u,u) = uuQ(u,u) , (150) 

where 

r(m)r(w) r(-w-ai)r(-m-q 2 ) 

^ ' j ~ r(i - iu) r(i - i«) r(i + »u + ai)r(i + i5 + a 2 ) ' 1 j 

This expression for Q(u,u) is symmetric to the substitution 

(u, ai,a 2 ) <-» (u,ai,a 2 ) (152) 
and can be written in the factorized form 

Q{u,u) — Q(u,ai,a 2 )Q(u,ai,a 2 ) , (153) 

where 

\ r(-iu- ai)r(-iu- a 2 ) 
Q(u,ai,a 2 ) = r 2 (l-m) (u) ' ( ] 



^^•j _ ""vi tu - ^177 °"H"i tu - "^7 (155) 



/ sin(7r(— iu — ai)) sin(7r(— iu — a 2 )) 
sin 2 (— i7ru 

The expression Q(u, ai,a 2 ) differs from the Baxter function in the holomorphic space 



Q ( u) = rHu rffr; u " 02) d56) 

only by the periodic function $(u) and therefore it can be considered also as a Baxter function. Note, 
however, that the function $(u) contains a square root singularity and, as a result, the recurrence 
relation for the function Q(u,u) differs from the similar relation for Q(u) by a sign in its right hand 
side 

(u — ia\)(u — ia 2 ) 
(u + i) 2 

To overcome this problem we can write Q(u,u) as follows 



Q(u + t,u) = -- " -Q(u,u). (157) 



Q(u, 5) =Q(«)Q(«) *(«.«). (158) 
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where the function <E> is given below 

= S in(.^ + S 1 ))sinKg + a 2 )) } 
sm(i7ra) smyinu) 

This additional factor <&(«., u) can be included in the definition of a new pseudo- vacuum state 

*o = ®(u,u)uu. (160) 

Really this pseudo-vacuum state can be considered as the additional factor for the wave function 
in the Baxter-Sklyanin representation providing correct hermicity properties of the hamiltonian and 
integrals of motion in this representation^ (see also Ref. [13]). We shall return to this problem in our 
future publications. 



8 Conclusion 

In this paper we established, that the gluon production amplitudes in the planar approximation could 
have the Mandelstam cut contributions in the multi-regge kinematics at some physical regions. For 
the cut corresponding to the composite states of n reggeized gluons the number of external particles 
should be k > 2 + 2n. The wave functions of these states in the adjoint representation satisfy the 
BFKL-like equation integrable in LLA and have the property of the holomorphic factorization. The 
corresponding holomorphic hamiltonian coincides with the local hamiltonian for an integrable open 
Heisenberg spin model. The Baxter equation for this model is reduced to a simple recurrent relation 
and can be solved in terms of the product of the T- functions. We constructed the wave functions of 
composite states of 2 and 3 gluons explicitly. 



I thank L. Faddeev, F. Smirnov, M. Staudacher, J. Bartels and A. Sabio Vera for helpful discussions. 

A Mandelstam cuts in planar diagrams 

Here we discuss an appearance of the Mandelstam cuts [33] in the crossing channels having adjoint 
representations of the color group SU(N C ) for the planar Feynman diagrams in the t'Hooft limit 
a <C 1, a c N c ~ 1 and calculate the impact factors corresponding to the multi-reggeon exchange. To 
begin with, let us consider the elastic amplitude A(s,t) for the gluon-gluon scattering in the Regge 
kinematics. It is well known, that in the leading logarithmic approximation the corresponding t- 
channel partial wave contains only one reggeized gluon pole. The contribution from the Pomeron 
exchange with color singlet quantum numbers is suppressed at large N c . The BFKL Pomeron appears 
as a composite state of two reggeized gluons and corresponds to the Mandelstam cut in the j-plane 
of the crossing channel. For the elastic amplitude the cuts in the adjoint representation appear in 
non-planar diagrams and are also suppressed at large N c . Indeed, according to S. Mandelstam these 
contributions should have the following form 

r d 2 k, r-sV'C-^) (- 5 y'(-(?-£) 2 ) 

^ s >*)~ h^t--^ 1 JZ rr 2 Mk±,q)Mk±,q), (A.l) 

J (2tt)^s fc2 (q-k) 2 

where j(t) are the Regge trajectories. The impact factors $ r are the integrals from the particle- 
reggeon scattering amplitudes f r (including the reggeon residues) over the invariants s r in the direct 
channel 

®r(k_L,q±) = / 7T~f r (Pr,k,q) , s 1 = (p A - k) 2 , s 2 = (pb + k) 2 . (A. 2) 

2 I thank Prof. F. Smirnov for discussions related to this important interpretation of the pseudo-vacuum state. 
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Here the integration contour L goes along the real axis above the right singularities of f r and below 
left ones according to the Feynman prescription. Only when the amplitude f r is constructed from 
the diagrams having both these singularities simultaneously the result of the integration is non-zero 
because in an opposite case we can shift the contour L from the real axis to infinity with a vanishing 
result. The Mandelstam cuts are absent also for the planar amplitude with five external particles. 

However, in the case of the six point amplitude there are planar diagrams in which the Mandelstam 
cuts are present. Let us denote the momenta of initial gluons by pa-,Pb and the momenta of final 
particles by pa> ,ki,k 2 , Pb' in an accordance with the order of multiplication of the corresponding 
color matrices T r . Then this cut appears in the physical region, where 

s = {pa+Pb) 2 > 0, si = {p A - +h) 2 < 0, s 2 = (h + k 2 ) 2 > 0, s 3 = {k 2 +p B ') 2 < 0. (A.3) 

This region corresponds to the transition of four particles with their momenta pa, —ki, —k 2 and 
pg to the two particles with the momenta pa< and pb 1 ■ In the multi-Regge kinematics, where the 
corresponding Sudakov parameters are strongly ordered 1 ^> — f3\ ^> — /3 2 , —on <C — a 2 <C 1, the 
integrands in the impact factors <& r 

$i(k,ki,q 2 ,) = / — — fi{p A ,k,ki,q 2 ) , $2(k,k 2 ,q 2 ,) = / -z— fz(PB, k, k 2 , q 2 ) (A.4) 

J L ^ m JL ^ m 

in the simplest case have only the poles in the integration variables a w 2kpA/s and j3 ~ 2/cps/s. 
Ill 1 



h 



(p A - k) 2 + it (fci +q 2 - k) 2 + ie - sa _ p + ie -sa0i - (h + q 2 - k) 2 + it 



1 1 _ 1 1 

' 2 ~ {PB + k) 2 + it (k 2 -q 2 + k) 2 +ie~ s p - k 2 + it sa 2 f3 - (k 2 - q 2 + k) 2 + it ' 

These poles are situated above and below the integration contours L due to the inequalities f3\ < 
0, a 2 < valid in the considered kinematical region where Si < 0, S3 < 0. Therefore the integrals are 
non-zero and can be calculated by residues 

^i(k,ki,q 2 ) = -= \ =-, $ 2 (k,k 2 ,q 2 ) = \ =- . (A.6) 

(fel + q 2 - k) 2 (k 2 ~q 2 + k) 2 

In the case of production of two gluons with the same helicity at the multi-Regge kinematics in 
the physical region where s\ < 0, s 2 > 0, S3 < the amplitude is proportional to the Born expression 

A 2 ^ = f 2 ^ 4 2s g T C A ) A -L g C(q 2 ,qi) T£ C1 9 C(q 3 , q 2 ) T* 2 g T£ B , (A.7) 
m\ w\ m\ 

where the Reggeon-Reggeon gluon vertex C is given above (see ([7])). In the lowest order approximation 
the corresponding proportionality factor f^o for the Mandelstam cut contribution in the £2-channel 
contains some additional multipliers from the effective vertices C in comparison with the above result 
(see gHH]) 



.g 2 N c /■ ^ d 2- 2efc ( _ S2)J (-P)-i ( _ S2)J (-(? 2 -fe) 2 )-i 

fL °= i — J — w w^w — 1 2 ' ( } 



where 



*i = + 92 - k) (q* 2 - k*) h,q 2 ) = k *l , (A.9) 

K[ + q 2 k 

$2 = qM - k* 2 - k*) (q 2 - k) $ 2 (fc, k 2 , q 2 ) = g3(g2 ~ \ . (A.10) 

q 2 -k 2 -k 
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In the weak coupling limit j = 1 and at e — » the amplitude f^o is 

lira / £G = nia fin ^ - i V a = ^ (4^) e , (A.ll) 

which coincides in this limit with the logarithm of the factor C introduced in Ref. [35]. This factor 
violates the Regge factorization of the BDS amplitude in the considered kinematical region due to the 
presence of the Mandelstam cut [38 . 

One can take into account the gluon reggeization in the channels t\ and using the following 
substitution in the above expressions 

' (sa - it) j ^- 2 , 1 r -» -(-s/3-ie) j ^- 2 . (A.12) 



— sa — k 2 + it s/3 — k 2 + it 

It would lead to the multiplication of the integrand with the real factor 

R= (-{k 1 + g 2 -k) \ f-fa-to+k) \ w ( _ ai fH)-i ( _ S3 )*3)-i . (A . 13) 

We can include also the diagrams with the reggeized gluon scattering in the crossing channel. It leads 
to the following expression for the Mandelstam contribution in LLA 

where G is the Green function satisfying the BFKL-like equation for the octet quantum numbers in 
i 2 -channel 

— — G(k, k', q 2 ; H-s 2 )) = K G(k, k', q 2 ; ln(-s 2 )) , 
o in s 2 

G(k, k', q 2 ; 0) = (2?r) „ <5 2 ~ 2e (fc - fc') . (A.15) 

Here the operator if in LLA can be expressed in terms of the Hamiltonian H which does not contain 
infrared divergencies 



H = 2 In lPl }yf +PiP* 2 ln| Pl2 | 2 -L + p\p 2 ln|p 12 | 2 J- , (A. 16) 



. . I" * +P*iP2 ln|pi 2 | 2 

i©r -P1P2 

where pi = k, p 2 = q — k. 



Let us consider now the Mandelstam cuts constructed from several reggeons. The non-vanishing 
contribution from the exchange of r + 1 reggeons appears in the planar diagrams only if the number 
of the external lines is n > 2r + 4. For the inelastic transition 2 — > 2 + 2r with the initial and final 
momenta pa-,Pb and pa> , fci, fe, fc 2r ,PB', respectively, the cut exists in the crossing channel with 
the momentum 

q = pa - pa> - = pb> - pb + = X! 9 ' ' ( A - 17 ) 

Z=l /=r+l i=l 

where q[ are momenta of reggeons forming the composite state. The corresponding amplitude has the 
form 

f r d 2 q' r ~'~ 1 ( sV^ — 

A(pA,PA', fel,..., k 2r ,PB>, PB) ~ / II^S 11- TM2 $1 ' " ^+l) $ 2 (<?! , •• ) • (A.18) 
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The impact factors $1,2 are given in terms of the integrals over the Sudakov parameters a[ — 
2q[pA/ ' s, 0[ = 2q[pB/ ' s from the reggeon-particle scattering amplitudes 

1=1 jl 

In QCD the tree expressions for fi_ 2 appearing in the planar diagrams are given below 
f =I _^ 1 1 1 

31 1 (pa <?i) 2 (pa - fco - <zi) 2 "" (pa - EU i - E[=o fcO 2 - E[=i 4 - E[=o ' 
/ 2 = / 2 11 



(pb + ?i) 2 (pb - W + q [r - (p fl + E[=i 4 - E?::+ 3 fcO 2 (p* + ELi 4 - h? ' 

where fco = pa 1 , fer+i = Pb 1 ■ The additional factors I\^ 2 contain effective reggeon vertices for the 
production and scattering of the gluons with the same helicity. They can be written in the multi-Regge 
kinematics JS]) as follows (cf. [I]) 

T _ t~t gt+i(Q ~ ELi V't ~ T! t Z\ k t) t~t a 

l\ (q* -Y!t\€ -y!^) Ir' 

where Q — Pa — Pas Q = Pb ~ Pb' and the Sudakov variables of the produced particles u\ — 
^iPa/s, ft — 2kipg/s are strongly ordered 

1 » |fti| » |/3 2 |... » Ifckl , |«i| « |a 2 | « ...|a 2fc | « 1 . (A.21) 
In these variables the functions /i i2 are given below 

h = h _A I I I ... 

— sa[ + it — s@ia[ - \Q — q[\ 2 + it —s0ia' 2 + it —s0 2 a' 2 — \Q - q[ — q 2 — ki | 2 + it 

f 2 = h _± 1 \ _J 

2 s0[ + it S a 2r p[ - \Q + qtf + ie sa 2r 0' 2 + * e sa 2 r-iP' 2 ~ \Q + q[ + q 2 - fori 2 +ie'"' 

where we took into account, that in the essential region of integration 

a t — , ft . (A.22) 

spi sa 2r -i+i 

In the physical region, where the signs of the Sudakov parameters of momenta fc; alternate with the 
index I 

01, a 2r < ; 2 , a 2r -i > ; 3 , a 2r _ 2 < ; ... , (A.23) 
which is equivalent to the following constraints on the invariants 

si < 0, s 2 < 0, s r < 0, s r+ i > 0, s r+2 < 0, s r+3 < 0, s 2r+ i < 0, s > , (A.24) 



the integrands in expressions for $12 contain poles above and below the integration contours L over 
y[, 0[. 



all variables a[, 0[. Therefore $i j2 are non-zero and can be calculated by taking residues from the 



poles in /tj 



$i(Sl,...,gL-i)=fT / , j q ' +1 p n ( A - 25 ) 

1=1 (Q* - El=i q' s * - El=i k t) (Q* - Etil q't* - Et=l K) 
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-^0 = II — ; %Wt^+i ; v'- 1 * i ' (A - 26) 

(=1 (Q + Ls=l Is - L s =l K2r-e.+l) (<9 + Lt=l 9t - Lt=l for-t+l) 

In the case of the production of 2r gluons with the same helicity the amplitude in N — 4 SUSY is 
proportional to the Born expression containing the effective Reggeon-Reggeon-gluon vertices C 0. 
The proportionality factor f^^2+2r for the Mandelstam cut constructed from r + 1 reggeized gluons 
can be written as follow 

In the leading logarithmic approximation the proportionality factor has the form 

f2— >2+2r 

Jlla 



rfl AT \ r r r ,,2e j2-2e„ ,,2e j2-2£„/ i 

^) / n ^i=^ n ^.^i *m-o»x *» n . ^ 

where we introduce the new notation pi for the reggeon momenta q[. The Green function satisfies the 
equation 



l-2e 



G(p,|/; Sr+1 )=](G(p,|/; Sr+ i), g(ft^;0) = n 2e ^""fo • (A.29) 

Here the operator _ftT in LLA can be expressed in terms of the Hamiltonian H which does not contain 
infrared divergencies 

K = w(t)-^H, w(t) = o0-ln^ , t = -\ q \\ 



g=m |pi1 y +x;g M+1 , (a.30) 
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where 



= In ]pi | 2 + In [Pi+il 2 + p« pf+i lnkl+i| 2 — ln|PM+i| 2 -^— • ( A - 31 ) 

P!P*+i PiP'+i 

Note, that the above hamiltonian has the property of the holomorphic separability 

r 

H = h + h*, h = ln^^+y hi,i +1 , (A.32) 

* (=1 

where 

/iu+i = In ^ + lnpi+i +p; ln/0;;+i h pi+i lnp M+ i . (A. 33) 

Pi Pi+i 

One can take into account also the enhanced contributions in the impact factors leading to the Regge- 
type dependence of the amplitude from other invariants Sj (i ^ r + 1). 



References 

[1] L. N. Lipatov, Sov. J. Nucl. Phys. 23 (1976) 338; 

V. S. Fadin, E. A. Kuraev, L. N. Lipatov, Phys. Lett. B 60 (1975) 50; 

E. A. Kuraev, L. N. Lipatov, V. S. Fadin, Sov. Phys. JETP 44 (1976) 443 ; 45 (1977) 199; 

I. I. Balitsky, L. N. Lipatov, Sov. J. Nucl. Phys. 28 (1978) 822. 



23 



[2] V. S. Fadin, R. Fiore, M. I. Kotsky, Phys. Lett. B 387 (1996) 593. 

[3] A. V. Kotikov, L. N. Lipatov, Nucl. Phys. B 582 (2000) 19. 

[4] L. N. Lipatov, Nucl. Phys. B 365 (1991) 614. 

[5] L. N. Lipatov, Phys. Lett. B 309 (1993) 394. 

[6] L. N. Lipatov, Sov. Phys. JETP 63 (1986) 904. 

[7] J. Bartels, Nucl. Phys. B 175 (1980) 365; 

J. Kwiecinski, M. Praszalowicz, Phys. Lett. B 94 (1980) 413. 

[8] L. N. Lipatov, Phys. Lett. B 251 (1990) 284. 

[9] L. N. Lipatov, Nucl. Phys. B 548 (1999) 328. 

[10] L. N. Lipatov High energy asymptotics of multi-colour QCD and exactly solvable lattice models, 
Padova preprint DFPD/93/TH/70, [hep^th/9311037[ unpublished. 

[11] L. N. Lipatov, JETP Lett. 59 (1994) 596; 

L. D. Faddeev, G. P. Korchemsky, Phys. Lett. B 342 (1995) 311. 

[12] H. J. de Vega, L. N. Lipatov, Phys. Rev. D64 (2001) 114019; 66 (2002) 074013-1. 

[13] S. E. Derkachov, G. P. Korchemsky, A. N. Manashov, Nucl. Phys. B617 (2001) 375; 

S. E. Derkachov, G. P. Korchemsky, J. Kotanski, A. N. Manashov, Nucl. Phys. B645 (2002) 237. 

[14] R. J. Baxter, Exactly solved Models in Statistical Mechanics, Academic Press, New York, 1982; 
E. K. Sklyanin, Lect. Notes in Phys. 226, Springer- Verlag, Berlin (1985); |arXiv:hep~th/9211111| 

[15] V. S. Fadin, L. N. Lipatov, Phys. Lett. B 429 (1998) 127; 
M. Ciafaloni and G. Camici, Phys. Lett. B 430 (1998) 349. 

[16] A. V. Kotikov, L. N. Lipatov, Nucl. Phys. B 661 (2003) 19. 

[17] J. M. Maldacena, Adv. Theor. Math. Phys. 2 (1998) 231. 

[18] S. S. Gubser, I. R. Klebanov, A. M. Polyakov, Phys. Lett. B 428 (1998) 105. 

[19] E. Wittcn, Adv. Theor. Math. Phys. 2 (1998) 253. 

[20] L. N. Lipatov, talk at "Perspectives in Hadronic Physics", Proc. of Conf. ICTP. Triest, Raly, 
May 1997. 

[21] V. M. Braun, S. E. Derkachev, A. N.Manashov, Phys. Rev. Lett. 81 (1998) 2020; 
A. V. Belitsky, Nucl. Phys. B558 (1999) 259. 

[22] A. V. Kotikov, L. N. Lipatov, V. N. Velizhanin, Phys. Lett. B 557 (2003) 114. 

[23] A. V. Kotikov, L. N. Lipatov, A. I. Onishchenko, V. N. Velizhanin, Phys. Lett. B 595 (2004) 
521; [Erratum-ibid. B 632 (2006) 754]. 1 

[24] S. Moch, J. A. M. Vermaseren, A. Vogt, Nucl. Phys. B 688 (2004) 101. 

[25] J. A. Minahan, K. Zarembo, JHEP 0303 (2003) 013. 

[26] V. A. Kazakov, A. Marshakov, J. A. Minahan, K. Zarembo, JHEP 0405 (2004) 024. 

[27] N. Beisert and M. Staudacher, Nucl. Phys. B 670 (2003) 439. 



24 



[28] A. V. Kotikov, L. N. Lipatov, A. Rej, M. Staudacher, V. N. Velizhanin, J. Stat. Mech. 0710, 
P10003 (2007). 

Z. Bajnok, R. Janik, T. Lukowski. larXiv:0811.4448l [hep-th], 
B. Eden, M. Staudacher, J. Stat. Mech. 0611 (2006) P014. 
N. Beisert, B. Eden, M. Staudacher, J. Stat. Mech. 0701 (2007) P021. 

Z. Bern, M. Czakon, L. J. Dixon, D. A. Kosowcr and V. A. Smirnov, Phys. Rev. D 75, 085010 
(2007) 

Z. Bern, L. J. Dixon, V. A. Smirnov, Phys. Rev. D 72, 085001 (2005). 
L. F. Alday and J. Maldacena, JHEP 0706, 064 (2007). 
L. F. Alday and J. Maldacena, JHEP 0711, 068 (2007). 

J. M. Drummond, J. Hcnn, V. A. Smirnov, E. Sokatchev, JHEP 0701 (2007) 064; 
J. M. Drummond, G. P. Korchemsky, E. Sokatchev, Nucl. Phys. B 795, 385 (2008). 

J. M. Drummond, J. Henn, G. P. Korchemsky, E. Sokatchev. farXiv:0709.2368l [hep-th] . 

J. Bartels, L. N. Lipatov, A. Sabio Vera. larXiv:0802.2 065 [hep-th]. 

J. Bartels, L. N. Lipatov, A. Sabio Vera. larXiv:0807.0894l [hep-th], 

Z. Bern, L. J. Dixon, D. A. Kosower, R. Roiban, M. Spradin, C. Vergu, A. Volovich, 
larXiv:0803. 14651 [hep-th]: 

J. M. Drummond, J. Henn, G. P. Korchemsky, E. Sokatchev, larXiv:0803.1466l [hep-th]. 

L. N. Lipatov, proceedings of the Conference "Quarks-2008", May 2008. larXiv:0810.3815l [hep-th], 

L. N. Lipatov, Scattering amplitudes in N~4 SUSY and open integrable spin chain, 
talk at the Workshop on Integrability in Gauge and String Theories, Utrecht, 10-15.08.08. 

[43] S. Mandelstam, Nuovo Cim. 30 (1963) 1148. 



25 



